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We propose a microscopic hamiltonian together with a master equation description to model 
stochastic hierarchical Random Telegraph Signal (RTS) or Pop-corn noise in nano-junctions. The 
microscopic model incorporates the crucial Coulomb correlations due to the trapped charges inside 
the junction or at the metal-oxide interface. The exact solution of the microscopic model is based 
on a generalization of the Nozieres-De Dominicis method devised to treat the problem of the edge 
singularity in the X-ray absorption and emission spectra of metals. In the master equation descrip- 
tion, the experimentally accessible transition rates are expressed in terms of the exact multi-channel 
Scattering matrix of the microscopic hamiltonian. 
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In the last two decades, the experimental as well 
as the theoretical investigations of noise in sub-micron 
structures have developed into a new and fascinating 
subfield of mesoscopic physics jy. Recently a resur- 
gence of interest in tunneling through quantum dots and 
other nano-structures has taken place, presumably due 
to progress in the field of low-dimensional strongly cor- 
related systems. As electronic devices are reduced to 
sub-micron sizes, stochastic fluctuations start to domi- 
nate their transport properties. In this letter we focus on 
generation-recombination noise that is often observed in 
semi-conductors and is interpreted as due to fluctuations 
in the number of carriers from the conduction (valence) 
band into traps in the band gap. Random Telegraph 
Signal noise (RTS) or Pop-corn noise is a special case 
of generation-recombination noise. The RTS may reflect 
the transition of single ions or trapped electrons between 
various meta-stable states inside the nano-structure. 

Tunneling through charge defects at interfaces or junc- 
tions has been discussed in the literature earlier |2| || . In 
particular, Matveev and Larkin showed that the current- 
voltage characteristic of a tunnel-junction with a local- 
ized impurity level exhibits a power-law singularity due 
to Coulomb interactions between the charge carriers in 
the leads and the local charge impurity H. Remark- 
ably enough, Geim et al. and Cobden and Muzykantskii 
measured the Fermi-edge singularity in resonant tunnel- 
ing and in RTS noise J|, |(|. More recently Doudin ct 
al. observed RTS fluctuations of up to 50% of the aver- 
age magneto-resistance in ultra-small junctions (tJ. Our 
investigations build further on these works. We do not 
discuss here the problem of dephasing as referred to in 

In this letter, we propose and solve a microscopic 
model as well as a phenomenological description of Pop- 
corn noise observed in tunneling through nano-junctions 
such as metal-oxide-metal nano-structures (extensions of 
this work to semi-conductor-oxide-semiconductor struc- 
tures, etc. are possible). At fixed chemical potentials, 



the current (conductance) jumps in a random yet hierar- 
chical manner between discrete levels (see Figure [l]). We 
argue that strong Coulomb correlations between trapped 
charges located inside the junction (or at the metal-oxide 
interface) and carriers in the leads govern the stochas- 
tic switching between the different values of the current 
(conductance). In the multi-channel and multi-trap case, 
we express the transition rates in a master equation de- 
scription in terms of the parameters of the microscopic 
hamiltonian. 

The model hamiltonian reads 
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where ct creates an electron of momentum k in lead 
a = L,R; more generally the index a may also label 
the various transverse channels a = l,...,n c inside the 
leads. For simplicity's sake we shall omit here spin de- 
grees of freedom, d- creates an electron in a localized 
state j inside the junction and at the interface. Traps 
can fluctuate in space between different mcta-stablc po- 
sitions or migrate, yet we shall neglect these fluctuations. 
e J ' 3 ' denotes the energy of and transfer between local- 
ized states. For simplicity, here we shall assume diagonal 
e?i = Sjj'Ej. e£ denotes the conduction band energies 
and \x a = e% F the chemical potential in channel a. More- 
over v a is the density of states at the a Fermi level. 
denotes the Coulomb repulsion between localized states 
inside the junction. Because of its crucial role in the hi- 
erarchical structure of tunneling, we shall treat this term 
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at the mean-field level ||. W kk , aa , represents the direct 
tunneling (non-diagonal matrix elements in the a indices) 
and backscattering (diagonal matrix elements) due to the 
metal-oxide interface. In the standard treatment of tun- 
neling, only the non-diagonal part of this term is taken 
into account. The physical effect of the backscattering 
potential is a renormalization of the momentum (energy)- 
dependent effective density of states in the usual formula 
for the tunneling current. W kk , aa , denotes the Coulomb 
scattering between localized and extended states and rep- 
resents the charge assisted tunneling. The terms W kk , aa , 
and W kk , aa , play a central role in this work and, should 
we emphasize, shall be treated exactly below. V ka de- 
notes the transfer integral between localized states inside 
the junction (or at the interface) and extended states in 
the leads. 

To give a simple physical picture of the tunneling with 
trapping we consider an hydrodynamical description of 
the electron fluids in the leads. At wave-lengths that 
are much larger than the average spacing between the 
particles, the charge and magnetization density fluctu- 
ations are described by harmonic oscillators that can 
be quantized. Coulomb correlations within the leads in 
the charge and spin degrees of freedom are included in 
this semi-classical heuristic argument. As charge SQ and 
magnetization 5M are being trapped in the junction (or 
at the interface), the electron fluids in the leads screen 
the accumulated charge and magnetization. After trap- 
ping has taken place, the energy increment consists in the 
charging and magnetizing terms proportional to SQ 2 , re- 
spectively, 5M 2 , plus a Coulomb term between the charge 
and spin fluctuations in the leads and the local charge 
and magnetization in the junction (or at the interface): 
the charge and spin density fluctuations can be described 
by displaced harmonic oscillators. Hence the overlap be- 
tween the de-trapped and the trapped configurations is a 
Gaussian in the charge and magnetization density fluctu- 
ations which vanishes in the infrared limit for Coulomb 
potential and one-dimensional longitudinal modes. In 
reality, the rate of charge SQ and magnetization SM ac- 
cumulation in the junction is proportional to the charge 
and spin current density differentials between the (chan- 
nels) right and left leads. As a result a dynamical cascade 
of Infrared Catastrophes occurs. 

There are a number of experiments |^|, ^, |^, [?J where 
the Markov property of the random switching process 
has been observed. Hence, we consider a description in 
real time in terms of a master equation. The idea is: 
i) to compute, for a fixed configuration a of trapped 
carriers, the charge assisted tunneling current {I) a at 
fixed chemical potentials ||. The configuration vector 
a = (ai,...,ajv) refers to cij = 1 if the local state 
j = l,...,N is occupied and dj = otherwise; ii) to 
evaluate the average current and the current-current cor- 
relation function through the junction by weighting each 



configuration appropriately || . Subsequently we give ex- 
pressions for the currents (/) a and the current-current 
correlation function. 

Combining a projection technique with partial tracing 
of bath (channels) degrees of freedom JTo| , we arrive at a 
master equation for the diagonal matrix elements of the 
local density operator, i.e., for the probability density at 
time t of being in a configuration a as 

p(a, t) = -A>(a, t) + J2 A >(a, i) , A a = ]T A? (2) 

a^a a^a 

with the transition rates from a to a given by 
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where p e ^ a represents the equilibrium density matrix for 
the bath (channels) with trap configuration a. Notice 
that a and a differ by a single site occupation. Further 
H' ah = P a H'Pb denotes a projected interaction picture 

of the transfer term H' = J2ka (^M C [-a4? + h.cA. The 
Markov approximation underlying equation (|^) requires 
the correlation time of the response function in equation 
(^) at finite temperature and in the presence of weak 
disorder to be much shorter than the characteristic time 
variations of the junction density matrix. For an arbi- 
trary number of localized states the master equation (|^) 
can be solved, in principle, by a method due to Kirchhoff 
pl| . A systematic Graph theoretic technique has also 
been developed by Weidlich |l2| . Figure [l] illustrates a 
three level hierarchical RTS noise. The crucial point here 
is that the transition rates, which are experimentally ac- 
cessible, can be expressed in terms of the local correlation 
function of the trapped states in the junction (or at the 
interface) . 

At this stage we have to tackle Dyson's equation for the 
transient propagator and from it infer the response func- 
tion (||). To solve Dyson's equation a substantial gener- 
alization of a method developed by Nozieres and De Do- 
minicis in the context of the X-ray edge problem in metals 
|i~3| , |l~4| is required. In the present case, Dyson's equation 
reduces to a system of singular integral equations that 
can be solved exactly by the techniques of Muskhelishvili 
and Vekua §, ||, (see also @). While the mathe- 
matics is rather involved let us attempt to explain in sim- 
ple terms the essence of the method. We consider the pro- 
cess of single electron transfer into and off the trap state 
jo by H', the occupation of all other localized states being 
fixed. The potential W* acting on the conduction elec- 
trons can be written as the sum of all W kk , (considered 
as matrix in the indices aa') for which aj = 1 plus the 
direct tunneling and backscattering term W kk , . A single 
electron transfer leads to Wg, where a differs from a by 
the occupation of the single site j - Immediately after 
this single electron transfer the new potential felt by the 
conduction electrons reads [Wl} kk , = [W^] kk , - [W a ] kk ,. 
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So the task can be decomposed into two steps, namely 
compute the propagator for a particle subject to the po- 
tential Wa, and then extract the transient propagator 
for the potential W~ which acts during a finite time 
interval. This transient propagator allows to obtain 
the response function required to compute the transi- 
tion rate (||). Following Nozieres and De Dominicis, we 
assume separable potentials W® k , aa , — Wq° Uk a Uk'a', 



an d V k a = V£ u ka where 
Uka is a cut-off function centered around the a Fermi sur- 
face, and take into account S-wave scattering only. This 
factorization simplifies the mathematics without affect- 
ing the physics, as it allows to perform readily the sums 
over k to obtain the local propagators at the junction. 
With these assumptions, the finite temperature transi- 
tion rate for an arbitrary number of traps, including di- 
rect tunneling, reads 



Al 
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2tt 
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cosh(/3(£f =F fi K )) + cos(we ai ) 
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where £o is of the order of the conduction bandwidths, 
T(x) denotes Euler's Gamma function and (3 = 1/fcsT 
the inverse temperature. The — sign refers to trapping 
(absorption) and the + sign to de-trapping (emission). 
k = 1, . . . , m c (< n c ) labels distinct chemical potentials 
p K . i?f = E a — E a with E a being the energy of the local- 
ized state renormalized (closed loop contributions) by the 
potential W a , i.e., E a — E a -\- AE° = J^f=i a 3 E 3 + AE°. 
The local Coulomb repulsion between trapped charges 
can be incorporated in the mean field approximation and 
contributes an additional term UN a (N a — l)/2 to the en- 
ergy E a , with N a the number of trapped charges in con- 
figuration a, and U the average potential in the junction. 
Furthermore A aiK is given by 



A aiK = -zJ^^A.+ti?!)- 1 )^ 



(BgDZlD i+ A% ia ,V$ (5) 



where jo denotes the trap subject to the change in occu- 



pation number (aj a ^ hj 



v a &aa' is the density 



of states at the Fermi level p a and we have defined the 
matrices (assuming summation over a) 



(AtfcW = [7r(tan6» Q ± i)u a ] 1 S aa > + W° 
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A* = [^r'x^W^K-i]- 1 ^)- 1 (7) 
(£* W = <W + ^(tan 9 a + ia(a, K))v a W£ a ' (8) 



with a(a, k) = +1 if p a < p K and —1 if p, a > p, K 
and (x K ) aa ' — S aa iS fiaif j iK . We furthermore have as- 
sumed that /ii > /J>2 > ■■■ > Hn c - The angle 9 a pa- 
rameterizes the short time behavior of the single par- 
ticle propagator in channel a Q. The matrix <Sf = 
D~+D a+ D^D a _ is related to the generalized multi- 
channel Scattering matrices, S a = D a+ D~^ and is di- 
agonalized by Bf, Sf = Bfe^^-Bf)- 1 , with the real 
eigenvalues <5f = diag((<5f )i, . . . , (<5f )„J. We use the no- 
tation e ai = ±2(«|) ai /7r - Tr((5|) 2 )/7r 2 . 

In the low temperature limit ksT <C \E£ =F \i ai | the 
transition rates show the power-law divergence with the 
exponents e ai known from the X-ray problem [T^, p"4|. 
The first correction in T is proportional to (fc.BT/|.E£ =p 
/i K |) 2 . In the high temperature limit, fc^T 3> \E? =F p ai , 
A shows a power- law behavior in T, 
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Specializing to a single trap, a single channel, and in 
the absence of direct tunneling the finite temperature ex- 
pression (|]) reduces to the result of § (their Eq. (5)). 
The effect of finite temperature is a broadening of the 
edge singularity, as illustrated in Figure ^[ A similar 
broadening occurs in the case of weak disorder in weakly 
coupled electrodes: The exponential decay of the elec- 
tron propagators in the electrodes leads to an imaginary 
contribution to the singular denominator, \E^±p K +i~/ a \, 

where 7 Q = 7ri/ a ndis,a^dis,Q and n dis,a, U d i s , a are the im- 
purity density and the weak disorder potential in the 
channel a, respectively. 

For a fixed configuration a it is possible to adapt the 
method provided by Caroli et al. to the calculation 
of the tunneling current in the presence of the backscat- 
tering W a a (a ^ a'). The result for the simplest case of 
two channels a — L, R is given by 



(/>• = 



cL, W^ R p h L (to)pU^) - W* L P h R {u)pi{u) 
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V r {uj)V a {uj) 
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where e is the electron charge, P ^(uj), p%(ui) represent the 
density of holes and electrons in the electrode a respec- 
tively, and 

= (1 - Wt L gT^m ~ W? R gl a {u)) 

Wi R W RL gT^)9 r A a ^l (11) 

with g^ a the retarded and advanced free Green's func- 
tions in channel a. As an illustration for the current- 
current correlation function, Sji(t), we choose a system 
with a single trap, a = a = 0, 1. In the master equation 
approach Sn becomes 



S H (t) = e-( A ?+ A o)l*l 



A?Aj 



(A? + AJ) 2 



((J>1 - </)o) 2 . (12) 
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In summary, we have been able to evaluate, for an ar- 
bitrary number of trapped charges, the transition rates 
of the master equation description of RTS noise and the 
tunneling current in terms of microscopic model hamil- 
tonian parameters. We emphasize that this approach is 
not restricted to two-level systems which is often found in 
the literature [|], but allows for several fluctuators that 
can create a more complex, in particular an hierarchical, 
RTS noise. In future work, we intend to take into account 
the correlations in the leads (channels) as well as spin de- 
grees of freedom. The latter are required, for example, to 
model RTS noise in magnetic nano-wires due to charge 
and magnetization trapping inside the junction. 
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FIG. 1: Computation of a three-level RTS current fluctuations 
using Eqs. (@) and @. 
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FIG. 2: Typical transition rate behavior in function of the 
energy of the localized state for the case of a single trap 
(Eq. (^)). The different curves correspond to the tempera- 
tures T = 0,Ti = T 3 1/3 » 40K,T 2 = T 3 2/3 « 1500 K,T 3 = 



HR/kt 



60,000 K. 



